library(sf)
library(ggplot2)
library(dplyr)
library(tmap)
Mekansal analiz ya da Mekansal Veri Bilimi (Spatial Data Science) verilerdeki mekansallığı dikkate alarak analiz eden, yani mekanin bilgisini analizin nesnesi haline getiren yeni bir bilim alanı
Hesaplama temelli istatistik, ekonometrik ve coğrafi bilgi kümesi olarak tanımlanabilir
-Mekansal analiz için en temel gereksinimlerden biri shape file denilen geometri dosyasıdır. - Basit anlamda shape file geometrik bilgiyi ve o geometriye ilişkin coğrafi özellikleri içerir. -Bir shape file’daki geometrik bilgiler nokta, çizgi ya da poligon ile temsil edilebilirler. -Uygulamada basitçe bir ülkedeki şehirlerin her birinin birer poligon olarak ifade edilebilir.
# turkey.shp ==> Shape File
# turkey.shx ==> Autodesk AutoCAD tarafindan derlenen, sekil ve font tanimlamalarini iceren dosya
# turkey.dbf ==> data base file
# prj,
#.gda,
#.sbx, vb.
https://data.humdata.org/dataset/turkey-administrative-boundaries-levels-0-1-2
turkey <- st_read("C:/Users/firat/Desktop/turkey_administrativelevels0_1_2/adm1a/tur_polbnda_adm1a.shp")
## Reading layer `tur_polbnda_adm1a' from data source `C:\Users\firat\Desktop\turkey_administrativelevels0_1_2\adm1a\tur_polbnda_adm1a.shp' using driver `ESRI Shapefile'
## Simple feature collection with 81 features and 8 fields
## Geometry type: MULTIPOLYGON
## Dimension: XY
## Bounding box: xmin: 25.66851 ymin: 35.80842 xmax: 44.81793 ymax: 42.10479
## Geodetic CRS: WGS 84
str(turkey)
## Classes 'sf' and 'data.frame': 81 obs. of 9 variables:
## $ adm1_tr : chr "ADANA" "ADIYAMAN" "AFYONKARAHISAR" "AGRI" ...
## $ adm1_en : chr "ADANA" "ADIYAMAN" "AFYONKARAHISAR" "AGRI" ...
## $ adm1 : chr "TUR001" "TUR002" "TUR003" "TUR004" ...
## $ Shape_Leng: num 8.27 5.49 7.15 6.59 4.43 ...
## $ Shape_Area: num 1.407 0.751 1.45 1.163 0.79 ...
## $ adm0_en : chr "TURKEY" "TURKEY" "TURKEY" "TURKEY" ...
## $ adm0_tr : chr "TÜRKIYE" "TÜRKIYE" "TÜRKIYE" "TÜRKIYE" ...
## $ adm0 : chr "TUR" "TUR" "TUR" "TUR" ...
## $ geometry :sfc_MULTIPOLYGON of length 81; first list element: List of 1
## ..$ :List of 1
## .. ..$ : num [1:2863, 1:2] 36.3 36.3 36.3 36.3 36.3 ...
## ..- attr(*, "class")= chr [1:3] "XY" "MULTIPOLYGON" "sfg"
## - attr(*, "sf_column")= chr "geometry"
## - attr(*, "agr")= Factor w/ 3 levels "constant","aggregate",..: NA NA NA NA NA NA NA NA
## ..- attr(*, "names")= chr [1:8] "adm1_tr" "adm1_en" "adm1" "Shape_Leng" ...
plot(turkey)
plot(st_geometry(turkey))
st_crs(turkey)
## Coordinate Reference System:
## User input: WGS 84
## wkt:
## GEOGCRS["WGS 84",
## DATUM["World Geodetic System 1984",
## ELLIPSOID["WGS 84",6378137,298.257223563,
## LENGTHUNIT["metre",1]]],
## PRIMEM["Greenwich",0,
## ANGLEUNIT["degree",0.0174532925199433]],
## CS[ellipsoidal,2],
## AXIS["latitude",north,
## ORDER[1],
## ANGLEUNIT["degree",0.0174532925199433]],
## AXIS["longitude",east,
## ORDER[2],
## ANGLEUNIT["degree",0.0174532925199433]],
## ID["EPSG",4326]]
Burada son satırda gördüğümüz “ID["EPSG",4326]]” bilgisi haritamızın projekte edilmiş olduğunu gösteriyor. Peki “st_crs” kodunu çalıştırdıktan sonra şöyle bir ifade ile karşılaşmış olsaydık ne yapacaktık?
#proj4string: "+proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs
Bu ifade bize haritamızın projekte edilmediğini gösterecekti çünkü “proj4string” ifadesini “+proj=longlat” takip ediyor. Eğer karşılaştığımız sonuç bu şekilde olsaydı, o zaman haritamızı projekte etmemiz gerekecekti. Peki hangi projeksiyon bilgisini kullanmamız lazım? Bu sorunun cevabi derin felsefi tartışmalar gerektirir ki bizim amacımız sadece haritamızı doğru bir şekilde görselleştirmek. Bunun için ülkelerin EPSG kodlarını bulmamız isimizi görecektir. Bunu yapabilmek için de ülkelerin UTM paftalarını bilmemiz gerekiyor.
Kısa bir Google araştırmasından sonra Türkiye’nin 4 farklı UTM paftasına sahip olduğunu görüyoruz.